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On  the  Theory  of  Flood  Waves  In  Rivers 

George  K.  Morikawa 
1.   Introduction 

Flood  waves,  in  rivers  vjith  small  bottom  slope  such 
that  subcrltical  flow  only  Is  encountered,  are  known  to 
progress  downstream  with  a  velocity  substantially  less 
than  the  propagation  velocity  of  small  amplitude  gravity 
waves,  (i.e.  the  downstream  characteristic  velocity). 
Such  is  the  case,  for  example,  in  the  Ohio  liiver  where 
a  goodly  amount  of  experimental  data  during  flood  condi- 
tions have  been  obtained  over  a  period  of  many  yeara  by 
Army  engineers.   On  the  basis  of  such  data  hydraulic 
engineers  have  devised  a  relatively  satisfactory  engineer- 
ing procedure  for  flood  prediction  and  control  (called 
flood  routing);  a  description  is  given  by  Gilcrest[l]. 
However,  the  recent  advent  of  high  speed  computing 
techniques  and  equipment  has  stimulated  renewed  interest 
and  effort  in  Improving  present  procedures  by  directly 
attacking  the  basic  non-linear  differential  equations 
governing  flow  in  open  channels  (cf.  NYU  reports  by 
Stoker  [2]  and  Isaacson,  Stoker  and  Troesch  [3]),   The 
present  report  is  concerned  with  a  study  of  the  signi- 
ficance of  the  propagation  velocity  of  flood  waves  and 
its  relationship  to  the  propagation  velocities  which 
characterize  the  basic  differential  equation  system. 

The  basic  differential  equations  for  flow  in  open 
channels  are  well-established.   In  the  present  investi- 
gation we  restrict  ourselves  to  the  simplest  possible 
form  of  the  equations,  namely,  flow  in  a  uniform  channel 
over  a  bottom  of  constant  slope.   Then  the  flow  system 
is  completely  described  by  two  equations  of  continuity 
and  momentum  for  the  channel  depth  h  and  the  flow  velocity 
V  as  functions  of  the  time  t  and  the  streamwise  distance  x. 
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2. 
(la)    h^  +  (hv)^  =  0 

(lb)    v^  +  vv^  +  gh^  +  ^ gS  =  0 

where  g  is  the  acceleration  of  gravity,  K  is  the  turbu- 
lence resistance  coefficient  (Che'zy)  and  S  is  the 
constant  bottom  slope.   Subscripts  denote  differentiation. 
The  turbulence  resistance  term  Kv  /h  can  be  written  in 

2  "2 

an  alternate  form  gv  /  2     where  c  =  (,gh)   is  the  signal 

yc 

velocity  and  y  =  l/K  corresponds  to  a  viscosity  coefficient. 

Thus  large  K,  corresponding  to  high  turbulence  resistance, 

implies  small  y  corresponding  to  small  viscosity.   The 

dimensxonless  velocity  ratio  v/c  is  called  the  Froude 

number,  based  on  a  characteristic  length  which,  in  open 

channel  flow,  is  the  depth.   In  the  analogy  with  compressible 

gas  dynamics,  v/c  is  called  the  llach  nuinber.   Although 

this  analogy  exists  we  emphasize  that  from  a  physical 

viewpoint  the  Froude  number  and  Mach  number  are  distinct; 

the  former  is  the  ration  of  inertia  to  gravity  forces 

and  the  latter  is  the  ratio  of  a  velocity  of  flow  to 

the  signal  velocity  of  small  disturbances.   (lb)  can  be 

refjarded  as  an  ep proximate  momentum  equation  valid 

under  tne  shallow  water  approximation  and  obtainable 

by  a  formal  mathematical  procedure  from  Euler's  equations 

of  motion,  e.g.  as  described  by  Friedrichs  [!|]  and  Keller 

[5],  with  the  addition  of  a  correction  term  for  turbulence 

resistance.      However,    vje   prefer    to  reiiard    (lb)   as   directly 

derivable  from  physical  considerations,  e,g.  as  obtained 

by  Stoker  [2],  where  the  dependent  variables  (h,v)  are 

considered  to  be  average  values  of  channel  depth  and 

flov;  velocity.   Then  the  turbulence  resistance  term  can 

be  reasonably  justified  on  the  basis  of  Prandtl's 
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3. 

"mixing  length  theory"  of  turbulence  [6].   The  impli- 
cation of  this  viewpoint  is  reserved  for  the  discussion 
in  Section  3.   In  application  to  actual  rivers,  the  basic 
equations  corresponding  to  (1)  include  additional  terms 
due  to  variable  channel  cross-section,  bottom  slope  and 
overbank  flow.   Also  other  forms  of  the  turbulence 
resistance  term,  which  seems  to  depend  primarily  on  tpe 
shape  of  the  channel  cross-section,  are  used  by  hydraulic 
engineers.   We  regard  these  factors  as  correction  terms 
which  will  modify  but  not  substantially  chanr^e  our 
conclusions  reached  on  the  basis  of  our  study  of  (1). 

The  basic  differential  equations  (1)  are  of 
hyperbolic  type  similar  to  those  of  compressible  gas 
dynamics,  a  subject  which  has  been  of  considerable  interest 
to  applied  mathematicians,  physicists  and  engineers  in 
recent  years.   The  usual  method  of  treating  this  tvpe  of 
equations  is  by  the  method  of  characteristics,  e.g.  as 
shown  in  [?]  and,  for  explicit  solutions  of  initial  and 
boundary  value  problems  for  these  equations,  by  finite 
difference  methods,  e.g.  as  discussed  in  the  IIYU  reports 
[2]    and  [31.   Kowe'"'er  because  of  the  presence  of  the 
turbulence  resistance  term  in  (1)  vje  find  that,  in  addition 
to  the  two  characteristic  velocities  (v  +  c),  a  propagation 
velocity  U  is  defined  where  v  g  U  <  (v+c).   U  is  the 
^relocity  of  the  known  steady  progressing  wave  solution 
(cf.  Section  2)  vihich  is  the  basis  of  the  usual  flood 
routing  procedures  (cf.  [31).   Previous  investigators 
studying  the  theory  of  flood  waves  ha-^'^e  noted  the 
significance  of  U  in  their  analyses.   Deymie'  [8]  has 
studied"  the  linearized  form  (perturbation  on  the  constant 
steady  flow  solution  (h  ,v  ))  of  (1),  resulting  in  a 
form  of  the  "telegraph  equation"  for  both  the  perturbed 
channel  depth  and  flow  velocity.   He  considered  a  problem 

G,  B,  Whitham  has  independently  studied  this  lineari- 
zation by  a  somewhat  different  analysis  but  with 
substantially  similar  results  in  an  unpublished  report 
at  I\[YU. 


'i!;  -iJ-T 


A  ri ..  5 


J  J  U  €•  i.  J  *  ■     ■ '  V  f" 


with  radiation  type  boundary  conditions  appropriate  for 

the  study  of  flood  waves:  given  zero  initial  conditions 

for  time  t  =  0  and  prescribed  channel  depth  or  flow 

velocity  as  a  function  of  time  at  a  fixed  point  x  =  0  for 

t  >  0,  to  determine  the  behavior  of  ther.e  perturbed 

quantities  downstream  (x  >  0)  and  for  t  >  0.   His  analysis 

shov;s    that    the    initial   disturbance   travels   dovjnstreara 

with  the  linearized  characteristic  velocity  (v   +  c  ) 

■'  ^  o    o 

but  that  the  main  or  maximum  disturbance,  which  initially 
travels  v/ith  this  velocity,  moves  downstream  v;ith  a 
velocity  which  asymptotically  (large  t)  approaches  the 
propagation  velocity  U  of  the  steady  progressing  wave. 
All  disturbances  are  completely  damped  eventually. 
Although  this  linearized  theory  yields  valuable  information 
on  the  behavior  of  certain  types  of  disturbance  wavea, 
it  is  not  completely  satisfactory  in  the  sense  that 
solutions  of  the  type  of  the  steady  progressing  wave 
are  not  included.   Recently  Hayashi  [9]  has  presented 
an  approximate  solution  of  the  initial  value  problem 
for  (1)  by  an  iteration  procedure,  but  it  has  been 
pointed  out  by  A.  Troesch"  that  this  solution  leads  to 
"breaking",  that  is,  bores  (or  shocks)  eventually 
develop,  even  for  disturbances  like  the  steady  progressing 
vmve.   In  this  report  these  questions  are  clarified. 
The  practical  applicability  of  our  study  to  actual  flood 
waves  will  be  included  in  a  future  report  by  A.  Troesch, 


Private  communication. 
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2 .   Approximate  Representation  of  Steady  Progressing 

VJave  and  Stretching  Tr ansf ormatl on . 

Since  one  of  our  primary  interests  is  in  obtaining 
approximate  time-dependent  solutions  (waves)  which 
travel  with  the  velocity  of  the  steady  progressing  wave", 
and  which  will  Include  the  SPVJ  solution,  we  transform 
(1)  to  a  coordinate  system  moving  downstream  with  this 
velocity  U,  namely,  (x,t)  — >(o,t)  where  o  =  x  -  Ut. 
Then  in  the  new  coordinate  system,  (1)  becomes 

(2a)    h^  +  [h(v-U)]^,  =  0 

Kv2 
(2b)    v^  +  (v-U)v^  +  gh^  +  £tL-  -   gs  =  0 

Setting  :r^  =  0  in  (2)  we  obtain  the  system  of  ordinary 
ot 

differential  equations  satisfied  by  the  SPW  solution 

denoted  by  h(6)  and  v(<3').   Prom  (2a)  we  determine 

h  (v  -U) 

(3a)     (V  -  U)  =  -2_° 

h 

which  we  substitute  into  (2b)  yielding  a  first  order 

ordinary  differential  equgtion  for  the  wave  profile 

H(c5-) 

(3b)    [c2(h/h^)^  -  v^d  -  U/v^)2]H^  + 

Kv^  |[1  +  (h/h^  -  l)U/v^]2  -  (h/h^)^/   =  0 

Hereafter,  "steady  progressing  wave"  is  denoted  by 
"SPW" and  refers  to  the  uniquely  determined,  bounded, 
non-trivial  solution  of  the  differential  equation  (3). 
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where  c  =  (gh  )  ,  h(oo  )  =  h  and  v(co  )  =  v  .   The 
propagation  velocity  U  is  found  directly  from  (2a)  for 

3 

(h  /h  )2  -  1 

ik)  U/v   =  --^^-^ 

(h,/h^  -  1) 

where  h(-cx)  )  =  h,  >  h  and 


(5) 


2     2 
1   _  '1 


Fig.  1 


The  shape  of  the  SPW  profile  is  shown  in  Fig.  1,   We 
note  here  that  the  constant  solution  h   satisfies  (3) 
and  the  constant  solution  v  satisfies  the  corresponding 

differential  equation  for  v(cr),  that  is,  i^Qt^^Q^i    i^  the 

steady  flov;  solution  which  satisfies  both  (1)  and  (2) 
with  the  condition 


(6) 


Kv. 


-  gS  =  0 


Both  the  steady  flow  and  SPV/  solutions  are  discussed  in 
greater  detail  in  Appendix  II  of  [3]. 

In  looking  for  approximate  time-dependent  solutions 
related  to  the  SPV  solution  an  ob'"^ious  generalization 
of  the  ordinary  linearization  of  Deymie'  [8]  is  to  carry 
out  a  perturbation  expansion  with  respect  to  a  snail 
parameter  e,  say,  on  the  known  steady  progressing  wa^e 
solution.   Putting 


fi 


oa^ 
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7. 
h((f,t)  =  h{(f)   +   eh^-'-'lo-.t)  +  e^h^^^(cJ,t)  +  ... 

v(d,t)  =  v(o)  +  ev(^^(c5',t)  +  e^v(2^(o',t)  +  ... 

in  (2)  and  equating  the  coefficients  of  the  first  power 
in  e,  we  obtain  the  first  order  system 

(7a)    hf^^  +  [(V  -  U)h^^^  +  hv^^h^  =  0 
(7b)   v(l'  +  ((¥  -  u)v<l'  +  gh'^')^  + 

h     V      h 

which  is  a  linear  system  for  (h   ,v^  ')  with  non-* 
constant  coefficients  depending  on  the  known  solution 
(h,v).   The  particular  case  (h,v)  =  (h  ,v  )  is  Deymie'' s 

first  order  linearized  system.   However,  explicit 
solutions  to  (7)  are  difficult  to  obtain  so  we  pursue 
another  approach. 

Returning  to  the  SPW  equations  (3),  ik)    and  (5), 
we  lool^  for  an  expansion  procedure  which  will  approximate 
the  SPW  solution  for  sriall  values  of  the  appropriate 
parameter  A  =  (h,  -  h  )/h  «1,  i.e.,  for  small  values 

of  the  waVe  height  (h,  -  h  )  with  respect  to  the  depth 

h  ,   So  we  let 
o 

(8)   h/h^  =  (1  +  h) 

where  h,/h  =  (1  +  X)  and  0  g  h  ^  X.   (i|.)  can  be 
expanded  as  a  function  of  X  only  yielding 

(9a)    U  =  U  +  XU-|_  +  \^U^   +  ... 


I  ■! 


i   I 
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or  specifically  for  the  steady  progressing  wave 

(9b)      u/v^  =la  .^.^^  ,..) 

Putting  (8)  and  (9)  in  (3)  and  retaining  only  "first 
order"  terms  in  h  and  X  gives 

(10)     h^  =  0 

which  yields  the  tri'^'^ial  constant  solution.   However, 
we  note  that  the  first  order  non-derivative  terms  are 
identically  zero  forcing  us  to  retain  the  second  order 
non-deri^^ative  terms.   Then  we  obtain 


(11a  )     h^  +  Ah(X  -  h)  =  0 


where 


(lib)     A  =i3Kv//Ul-i^)/{c^^   -  v^2/^) 


For  V  /c  <  2,  the  solution  of  (11)  is 
o'  o     '  ^ 


(12)     h(cJ)  =  1(1  -  tanh  ^) 

which  satisfies  the  required  properties  of  the  lowest 
order  approximate  solution  representing  the  steady 
progressing  v.'a'"'e  and  sEtows  that  the  deri'''ati'"'e  term 
h,  in  (11)  is  actually  of  second  order  in  X   and  that 

this  lowest  order  differential  equation  satisfied  by 
the  SPVJ  solution  is  non-linear.   (12)  shows  that  h  is 
a  function  of  (\o)  and  gives  us  the  motivation  of 
applying  a  stretching  transformation  '- '  =  Xo',  say, 
(cf.  Refs,  [i|]  and  [51)  with  an  expansion  procedure 


1  ; 


9. 

with  respect  to  the  parameter  X  in  order  to  obtain 
approximate  solutions  of  (2)  which  will  include  the 
SPVJ  solution.   It  is  now  clear  why  the  ordinary- 
linearization  procedure  of  Deyraie'  does  not  include 
the  3PW  solution.   In  the  following  two  sections  we 
shall  see  that  simultaneous  stretching  transformations 
on  both  6   and  t  yield  asymptotic  representations  of  the 
basic  differential  equation  system  (2)  and  consequently 
asymptotic  solutions  of  (2). 

2  -"• 
3 .   Asymptotic  Representation  for  (X,X  )   Stretching 

Transformation. 

A  priori  it  is  not  clear  from  the  basic  differential 
equation  system  (2)  what  type  of  stretching  on  the  time 
variable  will  yield  physically  reasonable  asymptotic 
representations.   It  appears  that  we  must  depend  rather 
heavily  on  our  knowledge  of  the  order  of  magnitude  of 
various  terms  in  (2)  based  on  physical  considerations, 
being  especially  mindful  of  the  appropriateness  of  (2) 
in  reasonably  describing  our  actual  physical  system. 
Thus,  the  stretching  transformation 


(13)    ^  =  Acf  ,  f=   X^t 


transforms  (2)  into 


(ll^.a)     Xh_   +  [h(v  -  U)l    =  0 


(li;b) 


xjxv    +  [(V  -  u)v    +  gh   A   +  H_  -  gs  =  0 
L  I  ^^  ^  J  h 


This  notation  is  used  to  describe  a  stretching 
transformation  on  {<S ,t)    — >(<U),7")  vjhere ''^-^  =  \(5   and 
r=  X^t. 


10. 

2 
Hence  the  ("^.  ,X  )  transformation  defined  by  (13)  implies 

that  in  (2)  the  time  deri^'ati^''e  terms  are  of  lower 

order  of  magnitude  (or  higher  order  in  \)    than  the 

space  deri^'^atl^e  terms  which  are  in  turn  of  lower 

order  of  magnitude  than  the  non-derivative  terms. 

Historically,  such  arfruments  on  orders  of  magnitude 

of  various  terms  in  the  Na^ier-otokes  equations  of 

motion  based  on  physical  considerations  led  to 

Prandtl's  concept  of  the  boundary  layer.   For  \  =  0 

in  Uk) 

=   0 


(15a) 

[h(v   -   U)] 

(l5b) 

■^f   -  63  =   0 

which  is  satisfied  by  the  steady  flow  solution 

(h,v)  =  (h  ,v  )   constant,  where  h  and  v   are  related 
oo  o      o 

by  (6).   Therefore  we  make  perturbation  expansions  with 
respect  to  the  parameter  X   of  (h,v)  on  the  constant 
solution  (h  ,v  )  by  expressing 

(16a)       h(d,t)  =  h^  +  nS^\c3  ,r)  +  \V^Uu),r)  +  ... 

(16b)    v(cy,t)  =  v^  +  \v(l^(tv,r)  +  X^v(2^(a),r)  +  ... 

Making  use  of  the  given  expansion  of  U/v  with  respect 

o 
to  \  we  put  (9a)  and  (16)  into  (lU)  and  equate 

coefficients  of  like-powers  of  A  in  the  usual  way. 

From  (I'^a)  we  obtain  for  consecutive  powers  of  X: 


'f : » , .'  .s 
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(17a)   X^:   [h^v^^^  +  (v^  -  U^)h^^h^  =  0 

(17b)  x':       [h,v(2).(v^-U^)h(2)]^^=.{h(}).(h(^)(v^l)-U,)]J 

(17c)   A^:   [h  v^^^  +  (v^-  U^)h(^^]    = 

and  so  on.   From  (li;b)  we  obtain  for  consecutive  powers 
of  X: 

(I8a)    X°:   Kv^^  -  gSh^  =  0 

(18b)    X^:   ^1 V-  =  0 

""o      ^o 

(18c)    x2:   h^(v^  -  Ujv^|)  ^  gh^h^^y 

(1)  2     p^(2)    .  (2) 
o         o       o 

(I8d)   X^:   h^(v„  -  U„)v^^2'  +  gh^hjj'  ^■ 

o 

_  ,^  ,  (1)  .  Vi  ,  (1)  ^  gh^^^h^^h 

and  so  on,   (l8a)  is  consistent  with  the  constant  steady- 
flow  condition  (6)  and  (17a)  and  (l8b)  implies  that 

U  /v  =  3/2  which  is  consistent  with  (9). 

°  °  (1)   (1) 

The  first  order  system  satisfied  by  (h^   ,v^   ) 

for  (X,X^)  stretching  is  obtained  by  (17a),  (17b), 

(l8a),  (l8b),  and  (l8c)  after  eliminating  the  second  order 


.     .   1   ■:> 
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combination  (Sv^^V'"  -  h^^Vh  )  ^  between  (17b)  and 

(I8c).   (I8b)  allows  a  separation  of  h^"^^  and  ^r^^^ 
resulting  in 

(19a,   .U>^.A{ch'l'u.hU'/VL-A4''.}  = 

where  A  is  given  by  (lib).   For  T<p— ,  =  0,  the  first 
order  SPW  solution  (12)  satisfies  (19a).   By  letting 

3v    h  (c^  -  V  2/1^) 

I19b)   V  -  -^^ ^^-^ 

o 

which  has  the  dimensions  of  kinematic  viscosity  and 
transforming  h^  '  by  a  simple  shift  factor, 

(19c)   u  =  I  v^(i  -  Sh^^Vh^) 
(19a)  becomes 

'"<"    ^%.j  -  "^"oj  ■  V  ° 

(19d)  is  a  non-linear  partial  differential  equation  of 
parabolic  type  (like  the  diffusion  or  heat  equation)  and 
has  received  much  attention  from  both  mathematicians 
and  engineers.   We  note  with  some  interest  that  \-ie   have 
obtained  the  identical  differential  equation  to  that 
obtained  by  previous  investigators  approximating  the 
Navler-Stokes  and  Reynold's  equations.   In  our  analysis 
we  have  started  directly  from  bulk  or  bydraulic-type 
considerations  of  continuity  and  momentum  with  a 
turbulence  resistance  term  having  some  basis  in  a  "mixing 
length  theory"  of  tui'bulence.   The  exact  significance  of 
this  observation  is,  of  course,  not  clear  since  so  many 


^y> :' 
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assumptions  ha'''e  produced  the  coincident  result  but  the 
connection  is,  at  least,  gratifying,   Kistorically  it 
appears  that  Bateman  [10]  first  studied  this  equation  in 
connection  with  the  problem  of  a  shock  wa"e  mo^''ing  in  a 
■viscous  compressible  medium.   Burgers  [11]  has  studied 
this  equation  in  his  well-known  contributions  to  the 
turbulence  problem.   Recently  Hopf  [12]  and  Cole  [13] 
ha^'e  independently  given  the  explicit  solution  to  the 
initial  value  problem  for  (19d)  and  -co  g  o  g  oo  , 
This  is  the  appropriate  problem  for  the  study  of  flood 
waves  and  the  desired  results  are  given  in  two  theorems 
by  Hopf,   The  first  theorem  states  that  if  u    (00)    is 
integrable  in  every  finite  Cj  -interval  and  that 

(20a)  60  "^  j  u^(4)dC  — >0  as  \cO\    — >  oo 


Then 

(20b)   u4^,  n 

where 


-00 


Jco  -^  ^^P  {  -  ^  Ffc..:,vi,r)j  an 

^    exp{-  -^  F(cO,^,r)]   d^ 


-00 


Vl 


2    r  ^ 
(20c)    F(6J, >/,?')  =  -^~^-  +  J  u^(4)d5 


is  a  regular  solution  of  (19d)  in  the  half  plane  f>  0 
that  satisfies  the  initial  condition 

U)  a 

(20d)   f  u(4,r)d4  -^  j   u^(4)d4   as  x  -^  a,  r-e>  0 

o  o 

for  every  a.   The  restrictions  on  the  initial  condition 
u  (/^))  at  V-   0  assures  both  uniqueness  and  existence 


i;c 


lk> 


of  the  solution  (19f )  for  all  T'>  0.   In  addition.  Cole 
has  considered  the  particular  problem  of  an  initial 
shock  condition""'  for  T"=  0: 


(20e)   u  ico)    =  -u,,  const,  for  x  <  0 


=  +Un 


for  X  >  0 


and  shows  that  the  desolations  from  the  SPW(12)  dampen 
out  like   exp(-u,  '^/l\.v).      This  initial  condition  and 
subsequent  3PV.|  at  'V  =   oo  are  shown  in  Fig.  2. 

h(l)t 


Fig.  2 

Hopf's  second  theorem  states  that  for  initial  data  in 
which  the  integral 


oo 


(21a)    M  = 


%(4)d^ 


-00 


exists  as  a  sum  of  two  improper  Integral  \ 


00       .  o 

O  -00 


'"  The  signs  for  x  >  0  are  reversed  from  Cole's  analysis 
to  correspond  to  our  physical  problem  of  flood  vm"es. 
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the  solution  ui/.J,?)    completely  dampens  out  for  large 

-  Vp  1/ 

time  like  V  alone  the  line  <''<-^/{2vT)  ' '^    =   const. 

For  example,  an  initial  condition  which  satisfies  (21a)  is 
(21b)    ^q(^0  =  +  u^,  const,  for  |x(  <  a 

=  0  for  |x|  >  a 

shov;n  in  Fig.  3. 

u(l) 


t  = 


n 


-a 


*=:^°M__ 


a 


Fig.  3 

ihe  behavior  of  the  solution  as  T^ — ^>  oo  for  this  "bump" 
type  initial  condition  is  consistent  with  the  behavior 
of  solutions  for  large  time  in  Deymie's  linearization. 

The  asymptotic  nature  of  the  stretched  and 
perturbed  system  (17)  and  (l8)  and  of  the  solutions  to 
the  corresponding  initial  value  nroblem  is  evident  by 
the  change  from  the  original  equation  (1)  and  (2)  of 
hyperbolic  type  to  appropriate  equations  of  parabolic 
type,   ^-'e  point  out  also  that,  in  this  approximation 
procedure,  the  parameter  X  has  a  definite  physical 
meaning  in  contrast  to  the  perturbation  parameter  of  the 
ordinary  linearization  (Deymie")  in  which  no  stretching 
is  involved;  and  that  it  is  essential  that  the  (X,X  ) 


a;  :'iii    . ' 
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stretchins  transformation  be  applied  in  the  moving  (cr,t) 
coordinate  system  whereas  the  (\,\)  stretching  trans- 
formation described  in  the  next  section  leads  to  the 
same  results  in  either  (x,t)  or  (Cjt)  coordinate  systems. 

The  second  order  system  satisfied  by  (h^   ,v^   ) 
for  (X,X  )  stretching  is  obtained  by  (17)  and  (l8)  after 
eliminating  the  third  order  combination  (2v'-^Vv  -h^-^Vh  )., 

between  (17c)  and  (l8d).   (l8c)  allows  a  separation  of 

(2)  (2) 

h^        and  v^  '  resulting  in 

-2h!i2)h(2)]   -  8  h(2)\  = 


.,.  l6Up    ,  (1)  2  ] 
+  A  ^[1)[ 2  ^  (£L_.)  ]y 


where  "l/v^  =  3/8  and  ^2/v^  =  -l/l6  for  the  SPW  (cf. 
(9)).   Unlike  U  ,  neither  U,  nor  Up  is  fixed  by  the 
approximation  procedure.   (22)  is  a  non-homogeneous 
linear  partial  differential  equation  of  parabolic  type 
with  non-constant  coefficients  as  functions  of  h^   (60,^) » 
given  by  the  solution  to  the  previous  first  order  problem. 

For  T^—   =  0  the  second  order  SPW  solution  of  the  form 

hr 

(23)   h(2)(^)  =  sech2  L^      fcosh^  ^^1     \   ^anh"  ^  ^p 

-^         n=o 

satisfies  (22).   We  vjill  not  carry  out  a  detailed  analysis 
of  the  initial  value  problem  for  (22)  in  this  report 
but  it  is  to  be  expected  that  the  solution  will  behave 
not  too  differently  from  that  of  the  ordinary  heat 
equation  since  the  dominating  terms  --  second  derivative 
term  with  respect  to  6J  and  first  derivative  term  with 
respect  to  time  r  of  h^ '   --  have  constant  coefficients. 
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[\. .   Asymptotic  Representation  for  (X,\)  Stretching 
Transformation. 

2 

The  approximation  procedure  with  the  (X,X  ) 

stretching  transformation  described  in  the  previous 
section  appears  to  be  the  most  natural  one  for  studying 
flood  waves  propaf^ating  downstream  with  the  velocity  of 
the  SPW,   However  we  shall  see  that  the  (X,A)  stretching 
transformation  studied  in  this  section  also  leads  to 

interesting  results:  1)  the  analysis  is  somewhat  simpler 

2 

than  that  v/ith  {\,\    )    stretching  -  a  consideration  which 

may  be  important  in  practical  application  to  flood  waves, 
2)  a  close  relationship  exists  between  (\,\)  and  (X,X  ) 
stretching,  and  3)  this  analysis  is  the  perturbation 
expansion  counterpart  to  Hayashl's  iteration  procedure 
and  has  the  advantage  of  being  consistent  with  the  SPi' 
solution. 

Thus,  we  consider  the  stretching  transformation 

(25)   ^O  =  XcT  ,  T"  =  \t 
which  transforms  (2)  into 


(26a)    h^.   +  [h(v  -  U)].    =  0 


(26b) 


x{v^  +  (V  -  U)v^  H-  gh^  j  +  111-  .  gs  =  0 


In  contrast  to  (X,A  )  stretching,  the  (A,X)  stretching 
transformation  defined  by  (13)  impliec  that  the  time  and 
space  derivative  terms  are  of  the  same  order  of  magnitude; 
but  like  (X,X^)  stretching,  (13)  still  Implies  that  all 
derivative  terms  are  of  lower  order  of  magnitude  as 
compared  to  non-derivative  terms.   As  in  the  previous 
section  we  make  perturbation  expansions  (16)  with  respect 
to  the  parameter  \   of  (h,v)  on  the  constant  steady  flow 


;•!.)  j 
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solution  (h  ,v  ).   Inserting  (19a)  and  (16)  into  (26) 
we  again  equate  coefficients  of  like-powers  of  X.   From 
(26a)  vie  obtain  for  consecutive  powers  of  \: 

(27a)  \^:    h[}^    +    [h  v^^^  +  (v   -  U  )h^^h   =  0 
C        o  o     O       ^J 

(27b)   \2.  ^{2)    ^    [h  v^^^  +  (v^  -  ^J^^^^l.  ■"    (h^^^^^^),,  =  0 
Y        o  o     o      i^  (j^) 

(27c)   X^:  h^^^  +  [h  v^^'  +  (v   -  U  )h^^h   + 
^  o  ^    o  o  Cj 

and   so  on.      From   (26b)    we   obtain  for   consecutive  powers 
Of   X: 

(28a)      X°:      Kv^^  -   gSh^  =  0 

(28b)      xS       {^- \-)    =  0 

o  o 

(28c)      X^:      v^i,)    .    (v^   -   Ujv^^l)    .   ghj/^)    ^ 

Kv„^        „(1)    2  ,„(2)         ,  (2) 


^    [(!_.)      ^    (2y h        )j    ^  Q 


^o  ^o  ^o  ^ 

(28d)      X^      v^^)    +    (v^   -   U^)v(j)    ^  gh^)    -f 

^[h(^)v(3).h^(v(^)..,)viJ). 
o  ^ 

O        /  ciV         V  .     '-V  n  \    _    ^^ 

"E"  ^ 2 "TT"      "^T     ~ 

O  V  0  0 


X  /  . 


f), 


I  :   .^^..^^ 


■      -1 

X 

T 


i:^':V   ^'i^- 


iA 
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and  so  on.   (28a)  is  consistent  with  the  constant  steady- 
flow  condition  (6). 

The  fir3t  order  system  satisfied  by  (h^-*- ^ ,  v^""- ^ ) 
for  (X,X)  stretching  is  obtained  by  (27a)  and  (28b) 
which  allows  a  separation  of  h^^  and  v^    resulting  in 

(29a)     h[i^    +  (3  ^    U^)hlV  =  0 
,',      do  o   uj 

If  the  propagation  velocity  U  is  taken  to  be  that  for 
the  SPW  then  U^/v  =  3/2  and  (29a)  becomes  simply 

(29b)    h^^^  =0 

(29b)  states  that  all  first  order  disturbances  remain 

unchanged  and  propagate  with  the  velocity  of  the  SPW, 

p 
We  note  that,  in  contrast  to  (X,\')  stretching,  (\,\) 

stretching  does  not  fix  the  value  of  U  although  the 

coefficient  (-^  v   -  U  )  of  hjj '  in  (29a)  implies  the 

importance  of  the  particular  value  U  /v  =  3/2.   Also, 
if  the  stretching  transformation  and  the  perturbation 
expansions  are  made  in  the  (x,t)  coordinate  system  as 
given  by  (1)  then  the  first  order  equation  for  (X,\) 
stretching  is 

(30)   h^i)  -i  V!^>^  =° 

which  is  identical  to  (29b)  under  the  coordinate  trans- 
formation (x,t)  — >(cr,  t)  where  cf  =  x  -  Ut  and 
U  =  u  =  3v  /2  is  the  first  order  approximation  to  the 

propagation  velocity  of  the  SPW,   Thus  the  following 
higher  order  systems  are  given  for  the  special  case  of  U 
given  by  (9b)  for  the  SPW, 


20. 


The  second  order  system  satisfied  by  (h^   ,v^  ') 
for  (A,X)  stretching  is  obtained  by  (27a),  (27b), 

jepj 


(28b)    and    (28c)    which  allows   a   separation  of  h^^    and 


(2) 
v'        resulting  in 

(30-  h^^    =   v[hiV    +  Ah^^^l    -   h^^Vh    )L  , 

where  A  is  defined  by  (lib)  and  v  by  (19b).   The  right 
hand  side  of  (30)  is  nullified  by  the  first  order  SPW 

(cf .  (19a)  for  -^^  =  0)  and  implies  h!i'  =  0  for  this 

special  case.   This  is  consistent  v;ith  the  original 
system  (2)  and  is,  in  this  sense,  an  improvement  over 
Hayashi's  iteration  procedure.   However,  any  first  order 

solution  which  does  not  nullify  the  rl£;ht  hand  side  of 

( 2 ) 

(30 )  implies  that  h    increases  linearly  with  time  and 

thus  becomes  unbounded  as  t  — > oo ,   This  implies  that 
the  (\,X)  stretching  transformation  inherently  leads  to 
solutions  singular  at  t  =  oo  for  almost  all  types  of 
disturbances  whether  a  perturbation  expansion  or  an 
iteration  (Hayashi)  is  made. 

The  third  order  system  satisfied  by  (h  '^  ,v^"^  ) 

for  (X,X)  stretching  is  obtained  by  (27)  and  (28)  and 

{  3  )      (  3  ) 
(28c)  allovjs  a  separation  of  h^    and  v    resulting  in 

(31)  h<3'  =  v(h'2UA(l-S^)h<2' 

o 


^^o.(l)  ^  3K,h(^^  2,    A   (l)ri    .h(^)  2.) 


CJ 


As  expected  from  the  study  of  the  second  order  system, 
i.e.,  equation  (30),  the  right  hand  side  of  (31)  is 
nullified  by  the  SPV,'  solution  (cf.  (23)  for  ^^^  =  0). 


21. 

This  result  is  again  consistent  with  the  original  system 

(2). 

If  we  novj  make  a  critical  comparison  of  the  results 

2 
of  our  analysis  by  (\,\)    and  (X,X  )  stretching,  our 

preference  from  the  viewpoint  of  the  theory  of  flood 

2 
waves  is  decidedly  in  favor  of  (X,\  )  stretching.   The 

2 
first  order  system  for  (\,A')  stretching  is  time- 
dependent  and  solutions  damp  out  as  t  — > oo  .   For  (\,\) 
stretching  we  must  investigate  the  second  order  system 
in  order  to  get  time-dependence  of  solutions  and  these 
solutions  are  unbounded  for  large  time. 
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Remarks 

In  our  study  of  the  theory  of  flood  waves  which 
propagate  downstream  with  the  velocity  of  the  steady 
propiresslng  wave  we  find  that  a  natural  way  of  obtaining 
appropriate  asymptotic  representations  is  by  stretching 
transformations  on  both  the  space  and  time  variables 
combined  with  the  usual  perturbation  expansion  but 
with  respect  to  the  stretching  parameter.   There  appears 
to  be  a  best  stretching  transformation  for  this  flood 
routing  problem,  namely,  the  (X,X  )  stretching,  which 
must  be  performed  in  the  appropriate  coordinate  system 
((cf,t)  instead  of  (x,t)).   However,  there  may  be  other 
stretching  transformations,  e.g.,  (\,A)  stretching,  which 
yield  additional  asymptotic  representations  giving 
valuable  information  and  which  may  be  particularly 
useful  in  application  to  other  physical  problems. 

Finally,  we  note  that  this  investigation  has  been 
made  for  the  particular  purpose  of  obtaining  reasonable 
asymptotic  representations  of  relevant  physical  phenomena. 
We  have  used  the  term  "asymptotic"  without  supplying  the 
estimates  that  are  implied  in  the  mathematical  usage  of 
the  word  "asymptotic".   ^  e  only  hope  further  work  will 
be  stimulated  in  this  direction. 


,it^ 
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